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'-^ I Abstract 

^ , The problem of critical behaviour of three dimensional random anisotropy magnets, 

which constitute a wide class of disordered magnets is considered. Previous results 
obtained in experiments, by Monte Carlo simulations and within different theoret- 
^ . ical approaches give evidence for a second order phase transition for anisotropic 

CN ' distributions of the local anisotropy axes, while for the case of isotropic distribution 

such transition is absent. This outcome is described by renormalization group in 
V^ . its field theoretical variant on the basis of the random anisotropy model. Consider- 

Cp I able attention is paid to the investigation of the effective critical behaviour which 

^-^ ' explains the observation of different behaviour in the same universality class. 
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rS ' 1 Introduction 



The influence of structural disorder on the critical behaviour remains to be 
one of the most attractive problems of phase transitions theory [1] . It is known 
that even small disarrangement in the structure of ideal physical systems may 
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have crucial consequences for their behaviour near a critical point [2,3]. In par- 
ticular, weak quenched disorder in magnetic systems can not only change the 
characteristics of the second order phase transition or the order of the phase 
transition into the low-temperature phase, it can modify also the nature of this 
phase, producing spin-glass order. Among disordered magnets depending on 
the type of randomness of their structure random site [2,3], random field [4], 
and random anisotropy [5] systems can be discriminated. All these magnets 
were extensively investigated, however, in our view, less attention was paid 
to the critical behaviour of the last class of systems i.e. random anisotropy 
magnets, in particular to random anisotropy effects on second order phase 
transitions. This motivates us to present an analysis of the problem by the 
field-theoretical renormalization group (RG) methods [6], which are standard 
for studying the second order phase transition of ideal systems; we also briefly 
review the present state of the investigations of the low-temperature phase of 
the random anisotropy magnets. The main goal of our study is to emphasize 
the crucial role played by distribution of the random axes in defining the uni- 
versity class of random anisotropy magnets, and hence the very possibility of a 
second order phase transition and its scenario. Moreover, in our study we will 
address the effective behaviour of random anisotropy magnets: it is frequently 
encountered in experimental and MC studies but less attention has been paid 
to it in the theoretical analysis. 

The paper is organized as follows: in the next Section we give a review of ex- 
perimental, numerical and theoretical analysis of the low-temperature phase in 
random anisotropy magnets. Since we apply in our study the field-theoretical 
RG, Section 3 discusses the functional representation of the model of a random 
anisotropy magnet for different types of distributions of the local anisotropy 
axes as well as the stability conditions for the effective Hamiltonians. Section 
4 presents the application of the field-theoretical RG approach to study the 
critical behaviour of the random anisotropy model (RAM). There, the RG 
functions are obtained in the two-loop approximations within the minimal 
subtraction scheme. The analysis of the asymptotic and the non- asymptotic 
critical behaviour is performed using the resummation techniques in space 
dimension d = 3. Section 5 concludes our study and summarizes the results. 



A brief review of earlier experimental, numerical, and theoretical 
studies 



Magnetic properties of random anisotropy magnets are described by the RAM 
[7]. This is a spin lattice model with each spin subjected to a local anisotropy 



of random orientation. The Hamiltonian reads: 

R,R' R 



Here, S'r is a classical m-component vector ("spin") on the site R of a rf- 
dimensional (hyper) cubic lattice, D > is the strength of the anisotropy, xr 
is a random unit vector pointing in the direction of the local anisotropy axis. 
The interaction Jr,r' is assumed to be ferromagnetic. For the short-range case 
it can be written in form of the nearest-neighbours interaction: 



J] 



R,R' 




|R, — R I — o; 
otherwise. 



where a is the distance between the nearest neighbours. The degree of disorder 
is controlled by the ratio D/J. Obviously, the random orientations are present 
in (1) only for m > 1. In the case m = 1 the random anisotropy term becomes 
constant and leads only to a shift of the Ising system free energy. 

The RAM (1) is relevant to describe a wide class of disordered magnets. It 
was first introduced to describe magnetic properties of amorphous alloys of 
rare-earth compounds with aspherical electron distributions and transition 
metals [7]. Today the majority of the amorphous alloys including rare-earth 
elements is recognized to be random anisotropy magnets. This class of amor- 
phous materials includes binary alloys [8,9,10,11,12,13,14], alloys of rare earths 
with two metallic components [15,16,17,18,19,20,21,22,23,24], alloys of rare 
earth with a metal and a non-metallic component [25,26,27,28,29] and mul- 
ticomponent alloys [30,31,32,33,34,35]. The structure of these compounds is 
characterized by an uniform (isotropic) distribution of the direction of the 
random anisotropy axes. Also crystalline compounds such as RExYi^x^h 
{RE = Tb, Dy) [36,37,38,39,40,41,42] are described by the RAM. However 
in this case, a distribution of the random anisotropy axes orientations is con- 
fined to several directions only. It turned out that random anisotropy plays 
a distinguished role also in the magnetic behaviour of certain uranium based 
amorphous alloys [43,44,45]. Although all materials mentioned above repre- 
sent the case of spin dimension m = 3, certain magnets are described by the 
RAM with m = 2 [46,47]. 

Beside traditional magnets, random anisotropy is present in such magnetic 
materials as molecular based magnets [48], nanocrystallity materials [49] and 
granular systems [50]. 

The problem of the nature of a low-temperature phase of random anisotropy 
magnets and peculiarities of the transition to this phase was investigated ex- 



tensively by theoretical approaches, numerical simulations and experiments. 
Below we summarize what has been obtained so far for the three dimensional 
systems. In a separate subsection, we mention also some results for the RAM 
at lower dimensions. 



2. 1 Experiments 



The earliest experimental data about amorphous alloys containing rare earth 
and transition metals are collected in the comprehensive review of Cochrane 
et al. [5] . Results of later experimental investigations are reviewed in the paper 
of Sellmyer and O'Shea [51]. Most of the investigated materials was recognized 
as systems with strong anisotropy. For many of them in the low-temperature 
phase an extremely large magnetic susceptibility was observed which was con- 
sistent with theoretically predicted peculiarities [52]. However no phase with a 
non-zero magnetization was obtained. The majority of investigations of amor- 
phous random anisotropy systems reports the low-temperature phase to be 
either a spin-glass phase or a so called correlated spin glass phase, in terms of 
the phenomenological theory [53] often used for the interpretation of data. 

Most of the examples given above concern systems with an isotropic distribu- 
tion of the orientation of the random anisotropy axes. However the crystalline 
compounds RE^Yi^^A^ (RE = Tb, Dy) [36,37,38,39,40,41,42] are examples 
of random anisotropy systems, where the distribution of local anisotropy axes 
is anisotropic. The rich phase diagram of these systems contrary to isotropic 
alloys presented evidence of a ferromagnetic phase for very weak anisotropy 
[38]. 

Therefore the experiments indicate the absence of a ferromagnetic phase for 
amorphous systems in general, with the possibility of ferromagnetism for crys- 
talline compounds with weak random anisotropy. The effective critical be- 
haviour of random anisotropy magnets was not considered. The question about 
effective critical exponents was touched only in Ref. [54], where the tempera- 
ture region of asymptotic critical regime was discussed for random anisotropy 
magnets. 



2.2 Monte Carlo simulations 



The majority of the computer simulations exploiting the RAM concerned 
the case of strong {D/J » 1) or more often infinitely strong anisotropy 
{D/J -^ oo). The earliest investigations report inconclusive results: both sta- 
bility [55,56] as well as instability [57,58] of the ferromagnetic order with 
respect to the spin-glass phase has been reported. However, data of later in- 



vestigations indicated the absence of ferromagnetism. The restriction to the 
infinitely strong anisotropy hmit led to a lack of long-range order in the ground 
state for m = 3 [59]. In this limit, the critical exponents at the transition 
to a low-temperature phase were found to be similar to those of the three- 
dimensional short-range interacting Ising spin glass [60] . The results of Monte 
Carlo calculations [61,62] confirmed an absence of ferromagnetism. 

On the other hand, Monte Carlo results for the RAM with m = 2 are consistent 
with the existence of a low-temperature phase of extremely large susceptibil- 
ity [61]. This feature was predicted theoretically for arbitrary m and weak 
anisotropy [52]. The possibility of the existence of such a phase was indicated 
also in investigations of spin models with p— fold random fields [62], which in 
the case m = 2, p = 2 correspond to the RAM with m = 2. 

Monte Carlo simulations for m = 2 show also the existence of a non-magnetic 
low-temperature phase with power-law decay of a pair correlation function 
predicted for the RAM [52]. The result [63] for infinite anisotropy is consistent 
with this quasi-long-range ordered phase. On the other hand such phase was 
not obtained for the model with p— fold fields dX J = D [64]. Moreover, another 
Monte Carlo study of the RAM with m = 2 and the same relation between J 
and D obtained critical exponents with values similar to the XF-ferromagnetic 
transition, except that the heat capacity critical exponent was found to be 
positive [68]. It was shown that the long-range order is not destroyed in the 
vortex- free model with p— fold fields [67] . 

The possibility of the existence of a quasi-long-ordered phase was also obtained 
for m = 3 in the case of weak anisotropy but assuming D/ J = oo for a part q 
of sites and D/J = for the rest of 1 — g sites [65]. Here, the quasi-long-range 
ordered phase appeared as an intermediate phase between the paramagnetic 
and the ferromagnetic ordering. However, for a non-zero anisotropy D the low 
temperature phase appears to be a quasi-long-range one [66]. 

To our knowledge, there exists only one the Monte Carlo study for the random 
axis anisotropic distribution. This is a study of a cubic model with random 
anisotropy, where the anisotropic axes are oriented along the edges of a cube 
[69]. For m = 2 the conventional XY second order phase transition to the 
ferromagnetic phase was found for weak random anisotropy, whereas a first- 
order transition to a domain type ferromagnetic phase was found for strong 
random anisotropy. For m = 3 both transitions were found to be of the first 
order. 

Therefore, although there is a certain inconsistency of Monte Carlo results, 
they mainly confirm the experimentally observed absence of ferromagnetism 
for random anisotropy magnets with isotropic distribution of random axes 
and bring about a second order phase transition for the weak randomness 



with anisotropic distribution. However, a low-temperature phase with power- 
law decay of correlation functions as observed in Monte Carlo simulations has 
not been found in the experiments so far. 



2.3 Theoretical treatment 



The first theoretical investigations of the RAM were performed within a mean 
field theory. Ferromagnetism was predicted [7,70] but the possibility of a spin- 
glass phase [71] was not excluded. Within the mean field approach the infinite- 
range interaction limit of the RAM was solved exactly [72] , indicating a second- 
order phase transition to ferromagnetic order. 

However, taking into account fluctuations lead the second order phase tran- 
sition disappeared. In the pioneer RG calculations performed for RAM with 
isotropic distribution of xr [73] no stable accessible fixed point of the RG trans- 
formation was found within the first order of £ = 4 — rf-expansion. Recently, 
this result was corroborated by two-loop [74,75,76] and even five-loop [77] cal- 
culations within the field-theoretical RG refined by resummation techniques. 
Meanwhile, an effective free energy was derived for large D. It was shown to 
have the same form as that of the random-bond Ising spin glass [78], demon- 
strating that the RAM can have a spin glass phase. Following the arguments 
of Imry and Ma [79] formulated for the random-field model it was shown that 
the d = 3 random anisotropy magnet should break into magnetic domains of 
size L ~ {J / DY [80] for weak anisotropy and thus no ferromagnetism was 
expected. Several different arguments were applied to the RAM [81] in order 
to demonstrate the absence of ferromagnetic order for space dimensions d < A. 
Although among these arguments the one for the limit m — *> oo [81] appeared 
to be erroneous [82], the lack of a ferromagnetism for the RAM with isotropic 
distribution of anisotropy axes for d < A was further supported by a Mermin- 
Wagner type proof using the replica trick [83]. The same result was obtained 
within the Migdal-Polyakov RG technique [84] avoiding the application of the 
replica trick. 

Investigating the equation of state of the RAM a zero magnetization and an 
infinite magnetic susceptibility were obtained in the low-temperature phase for 
any m [52]. Two-spin correlations in this phase possess a power law decay. As 
mentioned above, the phase with such a behaviour is called quasi-long-range 
ordered. The estimate of the susceptibility of low-temperature phase was cor- 
rected in another paper [85] and found to be finite with x ~ {D/J)~^ at rf = 3. 
A similar dependence of the susceptibility was obtained by other approaches 
[86,53]. The power law decay of spin correlations in the low-temperature phase 
was corroborated in particular for m = 2 [87]. But the last result is in dis- 
agreement with calculation results of Ref. [88], where spin-glass phase was 



obtained. Applying the functional renormalization group in 4 — e dimensions 
[89,90] to the RAM with m = 3, quasi-long-range order was found. 

Investigations of the RAM in the spherical model limit m — > oo were concen- 
trated on the question about the possibility of a spin glass phase. In this limit 
a ferromagnetic order was obtained for d > 2 and for D less than some critical 
value Dc, while for D larger Dc a spin glass phase was obtained for arbitrary 
d {Dc = for d < 2) [91]. For these investigations also the 1/m-expansion 
[86,92,93,94] was used. Applying the replica method a spin glass phase was 
found to exist below rf = 4 for arbitrary D [86]. But this spin glass solution 
later was shown to be unstable [92]. A stable spin glass phase for d < 4 was 
obtained avoiding the replica method [93]. However, the dynamics of spin glass 
order parameter for the RAM showed an instability of the spin glass phase 
[94] . These results found their confirmation in the mean field treatment of the 
m ^ oo limit [95], where the spin glass phase appeared only as a feature of 
this limit and no spin glass phase was obtained for finite m. 

The phenomenological theory [53] based on the continuous version of the 
Hamiltonian (1) and assuming correlations between randomly oriented an- 
isotropy axes turned out to be a more appropriate approach for the interpre- 
tation of the field dependence of the experimentally observed magnetization 
in the ordered phase. In this approach the spin correlation functions in differ- 
ent regimes of applied fields were analyzed [53] . In particular, the correlation 
length for small and zero fields was found to have the form C, ~ Ra{-^^y, 
where Ra is the correlation range for random axes. Such a phase was called a 
correlated spin-glass phase. 

The infinitely strong anisotropy limit of the RAM was investigated with the 
help of high-temperature expansions. Results of a Pade-analysis [96] indicated 
typical spin glass behaviour for space dimension d = 3, while in Ref. [97] it 
was concluded that the lower critical dimension for the RAM with m = 2 
is di = 3. The high-temperature series analysis of the Hamiltonian (1) on 
Cayley-trees [98] predicted ferromagnetic order, occurring for the number of 
nearest neighbours z > m and a spin glass order for z < m. 

All theoretical approaches summarized above concerned the RAM with iso- 
tropic random axes distribution. The anisotropic case was first investigated 
in a RG study [73] with a distribution of anisotropic axes, restricting direc- 
tions of the axes along the hypercube edges {cubic distribution). No accessible 
stable fixed point corresponding to a second order phase transition point was 
found. Investigating the RAM with mixed isotropic and cubic distributions it 
was found [99] that the presence of a random cubic anisotropy stabilizes the 
ferromagnetic phase. The possibility of a second order phase transition into 
a ferromagnetic phase with critical exponents of the diluted quenched Ising 
model for the RAM with a cubic distribution was pointed out in Ref. [100]. 



Such a behaviour was observed for a more general model [101] containing 
non-isotropic terms. It included the RAM with cubic distribution of random 
axes as a particular case. Subsequently, this result was corroborated within a 
two-loop RG calculations with resummation [75,76,102] done directly for the 
RAM and it was shown that the critical behaviour belongs to the universality 
class of the site-diluted Ising model. This result was further confirmed on the 
basis of a five- loop massive RG calculations [77]. 



2.4 RAM at low dimensions 



It is interesting to compare the above mentioned results with those for low- 
dimensional random anisotropy systems. For d=2 both theoretical approaches 
[103] and numerical simulations [104] lead to a spin-glass character of the 
low-temperature phase. An exception was a functional RG investigation [105] 
showing the instability of the ferromagnetic order and a logarithmically slow 
decrease of the correlation functions. 

Results for the RAM in d = 1 dimensions do not give such a consistent pic- 
ture. Investigations of the d = 1 infinity anisotropy limit [106] gave a spin-glass 
ground state. However, a Monte Carlo investigation for finite D [107] was un- 
able to find a ground state with zero magnetization. While another numerical 
study [108] demonstrated a decrease of the correlation length with D. Monte 
Carlo simulations [109] indicated a zero temperature magnetization. 



2.5 Some conclusions 



As one can see from the picture described above there exists a contradic- 
tory results concerning the transition of the random anisotropy magnets into 
the low-temperature phase. The majority of the results state that the d = 3 
RAM has no low-temperature phase with non-vanishing magnetization for 
the isotropic distribution of the local anisotropy axes. Investigations of the 
RAM with an anisotropic distribution of the anisotropy axes bring about the 
possibility of a ferromagnetic phase. Theory predicts a power-low decay of cor- 
relation functions in the low-temperature phase, which is confirmed by Monte 
Carlo simulations, however no phase with such features was observed in ex- 
periments. The question of effective critical behaviour of random anisotropy 
magnets remains unclear. 

Nowadays the field-theoretical RG approach serves as a reliable method to get 
accurate results describing critical behaviour [6]. For the RAM, only the papers 
[73,74,75,76,77,100,102] were devoted to such analysis. They concentrated on 
the study of the asymptotic critical behaviour of the RAM. Below we will 



consider results of these studies in more details, revisiting the c? = 3 RAM 
criticality within the minimal subtraction RG scheme. Doing so, we will study 
the effective critical behaviour of the random anisotropy magnets. It is this 
non-asymptotic critical behaviour that often is observed in experiments and 
Monte Carlo simulations of critical systems near second order phase transition. 
The analysis of the asymptotic and non- asymptotic critical behaviour is carried 
out in the Section 4, while in the next Section 3 the effective Hamiltonians for 
RAM with isotropic and cubic distributions of local axes are discussed. 



3 Effective Hamiltonians for random anisotropy systems 



The starting point for an RG analysis of the critical behaviour of the spin 
model (1) is an effective Hamiltonian, derived for a given random axes {xr,} 
distribution. The functional representation of the RAM and hence the effective 
Hamiltonian can be obtained using symmetry considerations [73,77], however 
it can be obtained also directly starting from the original spin Hamiltonian 
(1). Along the line of derivation, it will become clear, how the random axes 
distribution influences the symmetry of the effective Hamiltonian and hence 
defines the universality classes of our problem. The details of this procedure 
based on the Stratonovich- Hubbard transformation are given in the Appendix. 

The case we consider here corresponds to the non-equilibrium disorder [110]: 
variables ojr in (1) are randomly distributed and fixed (quenched) in a cer- 
tain configuration. As derived in the Appendix, the configuration-dependent 
partition function Z{{x-[i}) for the case of weak anisotropy D/J << 1 reads: 



Z{{x^}) ~ / (#)e^(^-'^^-), (2) 



where the integral means functional integration in the space of the m-compo- 
nent fields 0r = (0r . . . 0^" ) and the Hamiltonian 7Y(xr, 0r) reads: 



^oI^rI^ 



V0r' 



-Di (0r£r) +t;i|0R|'^ 



H(£R,0R) = -/rf'i?{^ 

z\$R\'{$iixjiy +...}. (3) 

Here, tq is proportional to the distance from the critical temperature and vi, 
Di, z are positive constants connected to the parameters of the spin Hamil- 
tonian (1) via relations explicitly given in the Appendix. 

This form of functional representation differs from the one obtained by symme- 
try consideration [73,77] by the presence of the term with coupling z, however. 



as it will be shown below, this term does not affect the critical behaviour of 
the RAM since the symmetry of the resulting effective Hamiltonian does not 
change when this term is omitted. In order to obtain the free energy F de- 
scribing the physical system for the quenched disorder [110] one has to average 
the logarithm of configuration dependent partition function Z{{x-r}) over all 
possible random configurations {x-r} of the directions of the anisotropy axes: 

-PF = J... J[d"'xn]P{{xn}) lnZ{{xn}). (4) 



Here, P({a;pt}) is the probability of a realization of a given configuration {xr}. 
If one assumes that there are no correlations between the directions of xr 
on different sites, the probability distribution is factorised into a product of 
distributions p{xu) of xr on each site R : 

^({£r}) = nP(^R)- (5) 

R 



To achieve the calculations of the logarithm of the partition function it is 
convenient to use the replica trick [83] introducing powers of Z{{xii}) which 
are easier to average: 



\nZ{{xR})=-^[Zi{xn})y 



(6) 

n=0 



Now the model must be completed by choosing a certain distribution p(xr) 
for the random variables xr. On the one hand, this distribution should be 
simple enough from the mathematical point of view, on the other hand, it 
must contain certain physical constraints. Aharony [73] considered two types 
of distributions of xr. The first is an isotropic one, where the random vec- 
tor xr points with equal probability in any direction in the ?72- dimensional 
hyperspace: 



Here r(a;) is Euler gamma- function, and the right-hand side presents the vol- 
ume of the m-dimensional hypersphere of unit radius. This distribution mimics 
an amorphous system without any preferred directions. The second distribu- 
tion restricts the vector xr to point with equal probability only along one of 
the 2m directions of the axes ki of a (hyper)cubic lattice: 



1 m 



2m ^^1 
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where S{y) are Kronecker's deltas. This distribution will be called cubic dis- 
tribution henceforth. The cubic distribution corresponds to a situation when 
an amorphous magnet still "remembers" the initial cubic lattice structure in 
spite of the random anisotropy or describes crystalline compounds with ran- 
dom but restricted anisotropy. It is an example of the cases of more general 
anisotropic distributions. 

Performing the average over the random variables {xr} for the isotropic dis- 
tribution (7) one ends up with the effective Hamiltonian [73]: 



n 



eff 



d'^Rl- 



^^lW? 



W 



+ uq\(p\ + 



Q,/3=l i,j 



"50f0f. 



(9) 



Q = l 



It describes in the replica limit, n ^ 0, the critical properties of model (1) with 
distribution (7). Here ^o = r^ — Di/m is the bare (non-renormalized) mass 



and Mo 



Df 



> 0, fo = f 1 + z/m > 0, Wq 



D? 



< are the bare 



m2(m+2) ^ "' ''U -i I -/"■' ^ ", "-U m(m+2) 

couplings. (f)f = ^f-^ are components of the n-replicated m-dimensional field, 
\0\^ = J2i J2a \(p?\ ■ The ratio of the couplings Wq and uq equals Wq/uq = —m 
and determines a region of physically allowed initial values in the {u — v — w)- 
space of couplings. 

For the cubic distribution (8) the average over the random variables {£r} 
leads to the effective Hamiltonian of the following form [73] : 



^e//=-/c? R<- 



2| ^|2 



V^ 



n 



+Moiv?r+^oE i<^"i'*+ 



m n r) m n 

«^oEE (C)'K) +1/0 EE 

i=la,/3=l ■ ' 



4 , „. y^ |j,a|4| 
a=l 

1=1 a=l 



(10) 



Again, one has to perform the replica limit n — > to get the physical quanti- 
ties. Here, the bare mass and couplings are determined as: ^iq = r^ — Di/m, 
uo = ^ > 0, vo = vi + z/m > 0, Wq = -1^ < 0. The last term in (10) 
has cubic symmetry. It does not result from the steps leading to "Heff but 
has to be included in (10) since it is generated by further application of the 
RG transformation. Therefore Hq can be of either sign. The symmetries of Wq 
terms differ in (9) and (10). Furthermore, although values of -Uq and Wq differ 
for Hamiltonians (9) and (10), their ratio again equals —m. 

It should be noted that for m = 1 the effective Hamiltonians (9) and (10) 
reduce to the traditional effective Hamiltonian of the Ising magnet with one 
coupling, demonstrating the degeneration of the random anisotropy term to a 
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constant for the Ising spins, as mentioned in the introduction. Indeed, in this 
case the Uq- and Wo-terms turn out to have the same symmetry (as well as the 
Vq- and i/o-terms in (10)). Moreover, due to the fixed ratio of these couplings 
they are equal in their absolute values but have opposite signs and therefore 
cancel one another. This leads finally to the usual 0*^-Hamiltonian. 

Other distributions leading to similar effective Hamiltonians as described abo- 
ve are discussed in the Ref. [77]. In particular, let the distribution p{x), have 
moments: Mi^i^^^i^^ = J d^x p{x) Xi^^Xi^ ■ ■ -^ik satisfying the conditions: 



M, 



ab 



m 



-Sab, 



M, 



abed 



A{5ab8cd + Sachd + 8adhc) + BSab^ac^ad-, (H) 



with Cauchy inequalities A{m + 2) + B > l/m and 3A + -B > l/rn^ for the 
parameters A and B. Then the effective Hamiltonian of the RAM with random 
anisotropy axis distribution with moments (11) turns out to have the form: 



n 



eff 



d^Rl- 



[ro-Z^i/m)|^|2+ V^ 



2711? 



a, 13=1 i,j 



BDi 



2 rn n 



EE 

i=la,/3=l 



-^oE 

a=l 



ia|4 



+ 



(12) 



However, a term YJlLi YZt=i (0?) has to be added here for the same reason, 
as it was included in (10). Thus the effective Hamiltonian corresponding to 
the distribution of J9(a;) with conditions (11) reads: 



^e//=-/c? R{- 



^^lw?+ 



V0 



n 



^oE E 

i=la,/3=l 



ka\'^ ( aI^ 



m n 

1/0 EE 



1=1 a=l 



Q = l 

n m 

n'+^oE EC0J0f0f, (13) 

o,/3=l i,j 



Such a Hamiltonian contains symmetry terms of both effective Hamiltonians 
(9) and (10). It was originally introduced in Ref. [100] to describe magnetic 
systems with single-ion anisotropy and non-magnetic impurities. The Hamil- 
tonian (13) was studied in the first order in e and no accessible and stable 
fixed point was found [100]. One should note that the i/q- term in the (12) 
as well as in (10) appears in a natural way, if one considers the original spin 
model (1) to contain a cubic symmetry term J2rJ2^i Sf^. 
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For distributions with A = in (11), the effective Hamiltonian (13) reduces 
simply to (10), since Wq = 0. As it was noted in Ref. [77] this feature is char- 
acteristic for all distributions, which orient the anisotropy axis only along the 
axes of an m-dimensional hypercube, including the cases when the anisotropy 
axis points in positive or in negative directions only. 

For the further analysis it is instructive to establish the regions in the space of 
the coupling constants where the free energy does not diverge in the ordered 
phase or in other words where the Hamiltonians (9) and (10) are stable. The 
regions vary for different ordered phases, described by certain nonvanishing 
values of the order parameter M. If one neglects field fluctuations and if the 
symmetry of the ordered phase is defined by (i) 0" = M then the regions of 
stability for the effective Hamiltonian (9) within the space of couplings are 
defined by the inequahty [100]: 



Vq + nuo + nwo > 0. (14) 

For the case (ii) 0" = MSaiSn the following condition has to be fulfilled [100]: 

vo + uo + Wo > 0. (15) 

The symmetries (iii) 0f = MSn and (iv) 0" = MSai lead to the same condi- 
tions (14) and (15) respectively. In a similar way four stability conditions are 
found [77] for the effective Hamiltonian (10): 

(i) mnuo + trvq + uwq + j/q > 0, (16) 

(ii) UQ + VQ + WQ + yQ> 0, (17) 

(iii) n-uo + "^o + ^""^o + ^o > 0, (18) 

(iv) muo + 'mvQ + WQ + yQ>Q. (19) 

As it was argued in Ref. [100], the legitimate stability conditions in the replica 
limit are those for which isotropy holds in the replica space. Then for Hamil- 
tonian (9) only the inequality (14) is a relevant stability condition, while for 
Hamiltonian (10) the appropriate conditions are inequalities (16) and (18). In 
the replica limit they reduce to conditions: 

vo > (20) 

and 

mvQ + I/O > 0, Vo + yo> 0. (21) 
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Together with the condition wo/uq = —m the above inequahties (20) and 
(21) determine all possible initial values of the effective Hamiltonians that 
can be used in the RG analysis, compared to some expiremental situation or 
computer simulations. 



4 Critical behaviour of the RAM as explained by the renormahza- 
tion group 



In this section the field-theoretical RG approach is applied to analyse critical 
properties of the RAM. Our previous field-theoretical RG studies exploited the 
massive renormalization scheme [74,75,76,102] and were restricted to investi- 
gation of asymptotic critical properties. Already in the two-loop approxima- 
tion [74,75,76,102] the critical behaviour was explained, then it was considered 
again in a five-loop study [77] confirming the former results. Here, we use the 
minimal subtraction RG scheme within the two-loop approximation to study 
the effective critical behaviour of the d = 3 RAM [111]. Since all essential 
features of asymptotic criticality become evident in this approximation, we 
consider it to be sufficient for analysing the effective critical behaviour. 



4-1 Renormalization 



The field-theoretical RG approach [6] is based on the renormalization of one- 
particle irreducible (HP) vertex functions defined as 



L N ^0 ^ ^ 

^(Ep^+E^.)ri^'''^(M,{n,/^o,Ko},Ao)=/e^^^^^"'+^^^'^^^^ 

i J 

< 02(ri) . . . cl)\rL)(p{Ri) . . . (P{Rn) >'^p(' dn... drLdRi . . . dR^, (22) 

where {ui^ = uo,vo,wo,yo, ...} is a set of bare couplings, /io is a bare mass, 
{p} and {k} are sets of the external momenta and Aq is the cut-off parameter. 
The angular brackets in (22) denote the statistical average over the Gibbs 
distribution with the effective Hamiltonian (9) or (10). 

The vertex functions appear to be divergent in the infrared limit Aq -^ oo. To 
remove the divergencies, a controlled rearrangement of the perturbative series 
for the vertex functions is performed. The finiteness of renormalized vertex 
functions T^^'^^ is ensured by imposing certain normalizing conditions. This 
leads to different renormalization schemes. In the following, the dimensional 
regularization with minimal subtraction is used [112]. 
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Normalizing conditions for the minimal subtraction scheme [112] are imposed 
at zero mass and have the following form: 



T('''\k,-k;fi;{u,}) 



k=0 



0, 



d 



r(°'2)(A;,-A;;/.;K}) 






j^2 — ij2 



fei-fc, = -^{4<S„-l) 



V^^'^\p■k,-k■fl■du^}) 



^fUi, 



1. 



(23) 



pk = -l/S/i^ 



^2 

Here, ui = /i^-^^-Uo are renormalized couplings and /i is the external momen- 
tum scale parameter, e = 4 — d. Expansions of the renormalizing factors of 
the field 0, Z^, the couplings {ui = u,v,w,y, ...}, Z^-, and for the operator 
0^, Zff,2, ensuring finiteness of two-point vertex function with one 0^ insertion, 
can be obtained from expressions (23). The functions T^:^^ are coefficients of 
terms of different symmetries in the expression for the four-point vertex func- 
tion r^^'"^)^^^^, containing the full tensorial structure. The last function for the 
isotropic distribution reads: 



p(0,4)ijfc« _-p{0,A) qijkl ,p(0,4)r.. p , p(0,4) /.ijfcZ 



'.P^Tl 



(24) 



while for the cubic distribution it is given by 



-p{0,A)ijkl _ p(0,4) qijkl _, p(0,4) C p i 



0/37T 



afSyr 



ijkl-'^ aP^T 



w 11 kl '-' ap'YT 1^ ■'' y -"■ iikl-'' amry 



(25) 



where 



^ ijkl ^ij^ik^ih 
1 



Sijkl — - i^ij^kl + Sik^jl + Sil^jk) , 
1 



^ijkJ — q {^ij^kl^aP^-iT+^ik^jl^a-i^PT+^il^jk^aT^p-y) 
"'ijkl^alSyr 

and S„h is the Kronecker's delta. 



Aaf3fT _ Q Q 
^ijkl — 'JiiklOn 



^ijkl 5 



(26) 



The renormalized vertex functions satisfy the following homogeneous RG equa- 
tion [61: 
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^(7<^^+70)}r(^'^)(W;{A;};/i;K}) = 0, (27) 

(except for the case L = 2 and A^ = 0) with the coefficients: 



Pu. = P^, (28) 

dhiZs , ^ 

d\nZ^2 
^^^ = -^h^' ^''^ 

determining the approach of a system to criticahty. The case L = 2 and A^ = 
fulfils an inhomogeneous equation and renormalizes additively. It is important 
for calculating the specific heat. Because of the scaling laws two exponent 
functions 7^, 7,^2 are enough to get all exponents. This also holds for the 
effective exponents above T^ in the approximation we use them here. 

The solution of the RG equation (27) via the method of characteristics intro- 
duces the flow equations: 

^?/ ■ 

£^ = PuA{n^})■ (31) 



The flow parameter i is related via a matching condition to the distance 
from the critical point r. In the limit £ — *> the scale-dependent values of 
the couplings Ui{i) will approach their stable fixed point (FP) values {u*}, if 
such a FP exists and if it is attainable from the initial values for (31) (initial 
"background" values of couplings) . 

The FPs of the system of equations (28) are defined as the set of real zeros of 
the /9-functions: 

/5«,(K})=0. (32) 



A FP is stable if all eigenvalues cjj of the stability matrix Bij = -q^ calculated 
at this FP have positive real parts. A stable FP corresponds to the critical 
point of the system only if this FP is reachable from physical meaningful ini- 
tial conditions. The critical exponents are then obtained from the 7-functions 
evaluated at this stable FP. The corresponding expressions for the exponent 
V of the correlation length and the exponent t] of the pair correlation function 
at Tc are expressed by the FP values of the 7-functions (29), (30): 
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^=[2-7^.(K})-70(K})]-\ (33) 

V = l4>i{<})- (34) 

The solution of the RG equation (27) together with the solution for vertex 
function r*^^'°^ proof scaling. Therefore, all other critical exponents may be 
obtained from the familiar scaling laws. For example, the susceptibility expo- 
nent 7 is given by: 

7=(2-77)z/, (35) 

and thus expressed by the 7-functions (29), (30) as: 

7-^1- ^*'";f» . (36) 

2-',4lu;}) 

Other exponents are found in a similar way. 
4-2 Asymptotic critical behaviour 



First, we analyse the asymptotic critical behaviour of the effective Hamiltoni- 
ans (9) and (10) within the renormalization scheme described above. For this 
purpose we write down the explicit expressions for the two-loop RG functions. 
The resulting FPs obtained by solving the system of the equations for the FP 
(32) are compared with the results of the massive scheme. This is the basic 
step before studying the effective critical behaviour in the next subsection 4.3. 



4.2.1 Isotropic distribution 

Imposing the renormalization conditions of the minimal subtraction scheme 
(23) to the vertex functions of the theory with the effective Hamiltonian (9) 
we found the [3- and 7-functions within the two-loop approximation (that is 
within the second and third order in couplings for the (3- and 7-functions 
correspondingly) : 

mn -|- 8 2 rn + 2 m+n+1 1 1 , 

Pu = —ue-\ u H uv-\ uw-\ — wv-\ — w — 

^ 6 3 3 3 2 

3mn+U . ll(m + 2) ^ 5 (m+2) . 

12 18 36 

(5mn+5m+5n + 87) , 11 (m-|-n-|-l) n 

^ -w'^u ^ -u^w- 

72 18 

5m-|-22 m+8 o^^ + ^ + Soln 

^g uwv —w V — w --V w, (37) 

17 



„ '"^ + 820 m + 5 

lJ^ = -ev^ — V +Zvu^ — vw- 

6 

3m+14 o llm+58 o 17m+5n+65 

V V u wvu — 

12 18 18 

5mn+82 o 13m + 56 o "^^ + 33m + 9n + 131 n 
^6—"" 1^^" 72 ^ "' (2^ 



m+n+4 2 o 2 

p^ = — ewH w +2'uwH — wv — 

6 3 

mn+3m+3n+17 o llm+lln+29 a 5mn + 82 o 

: w w u wu — 

24 18 36 

m+10 2 2 (m+5) 2 5m+34 ,^^, 

wv w V wvu, (39) 

36 9 18 ' ^ ' 



70 = + 



m + 2 m + n + 1 m + 2 

wv H uw H VU+ 

36 36 36 

mn+m+n+?> r, m+2 ^ mn+2 , , , 

7-^ ^i^ H — ;^^ H ^ — M^ (40) 

144 72 72 ^ ' 



2 '^nn + 2 m + 2 m + n+1 
^<^ 6 6 6 



w 



van + 22 ^Ti+2 m+2 2 ^rt + 2 

-u vu V vw— 

12 6 12 6 

m+n+1 ran + m + n + 3 2 

uw w . (41) 

6 24 ^ ' 

Here, e = 4 — c? and u, v, w are renormalized couplings. We recall that m is 
the spin component number and n is the number of replicas. The first line of 
each formula gives the one loop result in the corresponding RG function. We 
present the RG functions for any m and n, although in the further analysis 
of the RAM one considers the replica limit n = 0. However, the effective 
Hamiltonian (9) for non-zero n may have also other physical interpretations. 

An analysis of the RG functions (37)- (41) can be performed in two comple- 
mentary ways. First, one may apply the e-expansion [113]. Second, one may 
fix e (i.e. space dimension <£) in (37)-(41) and solve equations (37)-(39) for 
the FP numerically [114]. In first order of the e-expansion, no accessible sta- 
ble FP was obtained [73]. In second order, calculations performed within the 
massive scheme [74,75,76] did not change this picture qualitatively. Applica- 
tion of the ^-expansion within the minimal subtraction scheme leads to FPs 



coordinates which, being dependent on the renormahzing conditions [6] , differ 
from those obtained within the massive scheme [74]. The expansions of uni- 
versal quantities, e. g. critical exponents, coincide in both schemes. Leading 
to a quantitative information, these expansions for the disordered models [2,3] 
in general do not give reliable numerical values for the critical exponents. To 
give an example, in the case of a random Ising model the expansion turns out 
to be in y^, rather then e [115,116] and does not lead to reliable estimates 
at e = 1 (li = 3) [117]. Therefore we do not display here the results of the 
expansion in dimension e but we discuss in more details the results within the 
fixed dimension d approach. 

The weak-coupling expansions obtained within perturbation theory for the 
RG functions are asymptotic at best [6]. In order to obtain reliable results 
one has to apply resummation procedures. Here, the Pade-Borel resummation 
technique [118] is used. It consists in the following steps. For a given initial 
polynomial in several variables (in this case in three) for any finite series of 

p{u,v,w)= J2 ai,j,ku'v^w'' (42) 

l<i+j+k<3 



one introduces a "resolvent" polynomial in one auxiliary variable A by [119]: 
F{u,v,w;X)= Y. a^j^kuVw''y+^+''~\ (43) 

l<i+j+k<3 



with the obvious relation F{u, v,w; X = 1) = j3{u, v, w). Then, the Borel image 
of (43) is defined as: 

f«,,,„;.). j: ^-'f^^^-';:' . (44, 

i<i+j+k<3 \}+3 + f^ ly. 



This truncated series (44) is approximated by the Pade-approximant [1/1] (A) 
[120]. Then, the resummed /3-function is obtained from the formula: 

oo 

f^''^\u,v,w)= fdtexp{-t)[l/l]{t). (45) 



Applying the above procedure for the analysis of the RG functions (37)-(39) 
at fixed space dimension d = 3 and solving the FP equations (32) we get 
nine FPs. Table 1 lists the FPs for m = 2, 3, 4 in the region of interest: 
u* > 0, V* > 0, w* < 0. The numeration of FPs corresponds to that of 
Refs. [73,74,75,76]. FPs V and VII are not hsted in the Table 1 since their 
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FP 


m 


M* 


V* 


w* 


u 


7 


I 


\/m 











1/2 


1 




2 





1.1415 





0.684 


1.344 


II 


3 





1.0016 





0.720 


1.411 




4 





0.8878 





0.751 


1.470 


III 


Vm 


1.5281 








0.599 


1.185 




2 


0.1429 


0.9923 





0.665 


1.307 


IV 


3 


0.3411 


0.6965 





0.660 


1.298 




4 


0.4546 


0.5323 





0.656 


1.291 


VI 


2 


1.6637 





-1.8212 


0.452 


0.911 




3 


1.1351 





-0.9073 


0.470 


0.944 




2 


1.0126 


0.9058 


-0.6522 


0.631 


1.243 


VIII 


3 


1.0728 


0.7310 


-0.5483 


0.626 


1.235 




4 


1.1125 


0.6241 


-0.4823 


0.624 


1.230 



Table 1 

Fixed points and critical exponents u and 7 as obtained on the base of resuinined 
two-loop /?- and 7-functions for the RAM with isotropic distribution of random 
anisotropic axes. The only FP III is the stable one. 

coordinates are not in the region of interest. Performing the resummation 
procedure (42)-(45) sometimes one is lead to Pade-approximants containing a 
pole in the region of integration in (45). Then the principal value of the integral 
is taken. Coordinates of the fixed points obtained in this way are shown in 
the Table 1 by small numbers. Some of the FPs in Table 1 correspond to 
field-theoretical models, relevant for specific physical systems. In particular, 
FPs II-IV recover results for the isotropic m- vector, polymer {m = 0), and the 
diluted m-vector model correspondingly. All FPs in the Table 1 were already 
found within the e-expansion [73]. Their presence was further confirmed by 
the massive scheme of the field-theoretical RG [74,75,76]. Among these FPs 
only FP III is the stable one, however this FP is not accessible from the initial 
conditions for the coupling constants. 

Table 1 also presents the correlation length critical exponent u and the mag- 
netic susceptibility critical exponent 7 calculated for at the FPs. They are 
found from expressions u^^ = 2 — 7^ ({«*}) — 7^2 ({«*}) and 7~^ = 1 — 
7(^2({-u*})/(2 — 70 ({m*})), which are resummed in the same way as the /?- 
functions [121]. 

It should be noted that another stable FP which does not have an analogue 
among the FPs presented for the massive scheme [74,75,76] was found in the 
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region of interest. We do not list the FP in the Table 1 since we consider it 
uphysical for the reasons explained below. This new FP is not in correspon- 
dence with any FP found within e-expansion [73]. The presence of this FP 
can be an artifact of the resummation procedure caused by the nonlinearity 
of the system of the FP equations. Thus, the resummation can lead to the 
appearance of unphysical FPs. A possible check if a FP makes physical sense 
can be done by following the dependence of its coordinates on space dimen- 
sion d. If the non-trivial FP at d = 3 can be traced back to the non-trivial 
solution at (i = 4 (where the theory is known to contain only the Gaussian 
FP [6]) this may serve as strong evidence of the unphysical nature of such 
a FP. However the dimensional dependence of the coordinates of the FP in 
question could be calculated only up to d = 3.77, then no real solution can 
be obtained numerically. Moreover, starting from a certain value of rf, one has 
to take principal values of integrals in order to find the FP in the frames of 
the resummation procedure considered here. A similar problem was observed 
within two-loop approximation in the massive scheme [74] already at rf = 3. 
Furthermore, within the five loop order of the massive scheme [77] no evidence 
for this new stable FP was found. Therefore, we consider this new FP as an 
unphysical one. 



4-2.2 Cubic distribution 

Imposing the minimal subtraction renormalization conditions (23) to the ver- 
tex functions of the theory with the effective Hamiltonian (10) we get the 
following RG functions in the two-loop approximation: 



mn+8 ^ m + 2 n + 2 1 

3u=—su-\ u H vu-\ uw+yu-\ — wv— 

^ 6 3 3^3 

3mn + 14 o ll(m + 2) ,, 5 (m + 2) ^ 
12 18 36 

11 (n + 2) 2 5 (n + 2) o H 2 ^2 

u w uw u y y u— 

18 36 6 ^ 12^ 

3 5 -'- 2 -'- 2 ^ Ar\ 

-vuw wyu vuy v w w v, (46) 

2 6 6 9 9 ^ ' 



m + 8 2 o 2 

Py=—ev^ V +lvu^ — wv+yv— 

6 3 

3m + 14 o llm+58 2 5mn+82 2 ^+10 n 

V V U VU W V — 

12 18 36 36 

5 31 2 11 2 5 2 17 5n+34 
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Py^=—ew-\ w +2wu-\ — wv+yw— 

6 3 

3n + 14 o lln+58 ^ 5mn+82 2 

w w u wu — 

12 18 36 

31 2 m+lQ 2 11 2 5 2 
-w v-^^v W--W y--y w- 

17 5m + 34 5 

-—wuy — — wvu--vwy, (4J 

18 



3 4 

Py = -ey + -t/^ + 2yu + 2?/t; + 2wy + -wt;- 

17 , 5mn + 82 2 23 2 23 o 23 2 

— y u y y u y v y w— 

12^ 36 ^6 6 6 

5m + 82 2 5n + 82 2 ^+18 2 m + 18 2 

V y w y w V v w— 

36 36 9 9 

41 5n + 82 5m + 82 8 

—vwy — — wyu — — vuy--vuw, (49) 



7-^ = 



mn + 29 y^ n + 2 , m + 29 yu 

u^+—-\ w H v^+- — h 

72 24 72 72 12 

wv yv m + 2 n + 2 wy 

vu-\ wu-\ (50) 

12 12 36 36 12 ^ ^ 



2 iTin + 2 m + 2 n + 2 y 



-w -\ 

6 6 6 2 

mn + 2 ^ m + 2 ^ n + 2 ^ y 
-u — — V —^w 



2 



12 12 12 4 

m + 2 n + 2 yu wv yv wy ,__, 

vu wu . (51) 

6 6 2 2 2 2 ^ ^ 

Here, u, v, y, y stand for the renormalized couplings and again the rephca hmit 
n = will be considered. 

As noted in the previous subsection 4.2.1, two independent ways to analyse 
expressions (46)-(51) are possible. The e-expansion to first order in e was 
performed in Ref. [73]. It reported 14 FPs but no stable and accessible one 
was found among them. However, subsequently Mukamel and Grinstein [100] 
pointed to the possibility of a second order phase transition with exponents 
of the weakly diluted quenched Ising model [115]. Indeed, performing pertur- 
bation theory expansion to order e^ not only the corrections to the existing 
FP are obtained but new FPs [75,76,102] appear. Their appearance is caused 
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by the degeneracy of /3-functions at the one loop leveL In particular, in the 
replica limit the one- loop parts (the first lines respectively) of /9^, Py (48), 
(49) and /9„, Py (46), (49) are degenerated for any m at u = v = and at 
w = V = 0, correspondingly. Another one-loop degeneracy of the /?„, P^, Py 
functions with w = occurs aX m = 2. The last degeneracy was observed also 
for the diluted cubic model [122,123]. This implies a y^-expansion [115,116] 
for the FPs [75,76,102] rather than an 5-expansion. Among the FPs found 
with the help of the y^-expansion, the FP with coordinates w* < 0, y* > 0, 
■u* = t>* = is stable. Moreover, this point is reachable from the initial values 
of the couplings. As far as it is the FP of the diluted Ising model one con- 
cludes, that in the critical region the RAM with a cubic distribution of random 
anisotropy axes (10) decouples into m independent diluted Ising models and 
the phase transition is governed by the familiar random Ising model critical 
exponents. 

This picture found its further confirmation by the fixed d approach in the two- 
loop [75,76,102] and later in the five- loop [77] approximation of the massive 
scheme. The minimal subtraction scheme corroborates these results. Applying 
the resummation procedure (43)-(45) to the /^-functions (46)-(49) we get 16 
FPs. In Table 2 we present the numerical values of the FPs coordinates with 
u* > 0, V* > 0, w* < 0. The numeration of FPs correspond to that of Refs. 
[73,75,76,102]. Only FP XV and the "polymer" FP III are stable. However the 
last one is not reachable from physical initial conditions. FP XV in Table 2 
corresponds to the stable FP found in the y/e-expansion [102], and has coor- 
dinates with u* = V* = 0, w* < and y* > 0, thus it is accessible from typical 
initial values for the couplings. 

Applying the resummation procedure (43)-(45) to the /3-functions (46)-(49) 
we have not found any other stable FPs in the region of interest. Thus one 
concludes that the effective Hamiltonian (8) in the asymptotic region reduces 
to a product of m effective Hamiltonians of a weakly diluted quenched random 
site Ising model. This means that independent on the dimension m of the order 
parameter field in (3) the system is characterized by the same set of critical 
exponents which are those of the weakly diluted random site quenched Ising 
model. 

It is important to note that the FP XV does not describe critical behaviour 
if the random anisotropy axes distribution differs from the cubic one. In such 
a case the effective Hamiltonian (13) has to be considered and FP XV turns 
out to be unstable with respect to the term with zq. It was stated already 
in the one- loop study [100] and checked by nonperturbative arguments [77]. 
However critical properties of Hamiltonian (13) were studied only in one-loop 
order so far. 

From the stability exponents of the stable FP XV, one can obtain informa- 
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FP 
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II 
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1.411 
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0.8878 
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1.470 


III 
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1.5281 











0.599 


1.185 


V 


Vm 











1.3146 


0.644 


1.269 




2 


0.1429 


0.9923 








0.665 


1.307 


VI 


3 


0.3411 


0.6965 








0.660 


1.298 




4 


0.4546 


0.5323 








0.656 


1.291 


VII 


Vm 


2.5382 





-2.5382 





1/2 


1 




2 





1.9719 





-1.3146 


0.644 


1.269 


VIII 


3 





0.8568 





0.2270 


0.717 


1.405 




4 





0.4435 





0.7555 


0.715 


1.403 




2 


-0.2506 


2.4494 





-1.6330 


0.667 


1.312 


IX 


3 


0.3404 


0.7275 





-0.0511 


0.660 


1.298 




4 


0.4325 


0.4342 





0.2149 


0.657 


1.292 


X 


Vm 


0.7311 





-0.7311 


1.3146 


0.644 


1.269 


XV 


Vm 








-0.2506 


1.6330 


0.667 


1.312 



Table 2 

Fixed points coordinates and critical exponents v and 7 as obtained on the base 
of resummed two-loop (3- and 7-functions for the RAM with cubic distribution of 
random anisotropic axes. FP III and FP XV are stable. Note that the exponents 
may take the same values for different FPs. 

tion about how fast the couplings approach their asymptotic values when one 
goes to Tc- Approaching the critical point the correction to scaling enters the 
power-law behaviour of thermodynamical quantities. This correction is de- 
fined by the smallest stability exponent of the stable accessible FP. We get 
the following values for the stability exponents of the FP XV: uji ~ 0.1935, 
UJ2 ~ 0.4939, oju ~ 0.0036, oj^ ~ 0.0036. Here, uii and UJ2 are the familiar sta- 
bility exponents of the random-site Ising model. Currently, accurate values of 
these exponents within five-loop [117] and six-loop [124,125] approximations 
are known. Let us note, that the leading correction to scaling exponents for 
the model (10) is determined by the stability exponents a;„, uiy. General ar- 
guments [77] state that the leading correction exponent can be expressed by 
oj = oJu = ojy = —arjvr-, whcrc ar and v^. are the random-site Ising model 
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exponent of the specific heat and the correlation length respectively. Indeed, 
this is fulfilled: a.r/vr ~ —0.0036, in two loop order. The Wegner correction to 
scaling A = ujVr = ~oit is equal to the random-site Ising model specific heat 
critical exponent with opposite sign [77], which has been calculated with higher 
accuracy up to six loops [2,117,126] and is estimated as A ^ 0.049 ± 0.009 
[77]. 

For the other FPs, we recover the familiar two-loop numerical results for the 
Gaussian (FPs I, VII), the m- vector (FP II), the polymer 0(n = 0) (FP III), 
the Ising (FPs V, X), the diluted m-vector (FP VI), and cubic (FP VIII) 
models. FP IX belongs to a new universality class. In the Table 2 we give the 
numerical values of the critical exponents in these FPs as well: if the flow from 
the initial values of couplings pass near one of these FPs one may observe an 
effective critical behaviour governed by the values of the corresponding critical 
exponents. This will be discussed in the next subsection in more detail. 



^.c? Effective critical behaviour 



It is hard to approach the critical point close enough to observe asymptotic 
critical behaviour in the experiments as well as in the computer simulations. 
Therefore, very often one deals with the nonasymptotic critical behaviour de- 
scribed by the effective critical exponents. These exponents are defined as log- 
arithmic derivatives of corresponding thermodynamic quantities with respect 
to the reduced temperature r [127]. Within the RG approach the effective crit- 
ical exponents are defined in the region, where the scale dependent couplings 
Ui{i) have not reached their FP values and depend on the flow parameter 
i. In particular for the susceptibility critical exponent '-fcs it results in the 
expression: 



%VW = l- ,^'*'™' +.... (52) 

2-7<^(K(^)}) 

where the dots mean contributions from amplitude parts of the corresponding 
thermodynamic function. In our calculations we neglect these parts assuming 
contribution of the amplitude functions to the crossover to be small [128]. 

If the system of differential equations (31) has a stable FP {u*} accessible from 
the initial conditions, this FP corresponds to a critical point of the system. 
Then RG flows in the limit £ — > reach this FP, and the effective critical 
exponents in this limit coincide with their asymptotic counterparts. However 
the initial conditions being near the origin in the space of couplings (i. e. 
near the Gaussian FP) the flow may approach unstable FPs on its way to 
the stable one. Then the effective exponents may take values corresponding 
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to these unstable FPs. 

In the next subsections the effective critical behaviour of the both Hamil- 
tonians (9) and (10) is investigated in the two-loop approximations of the 
minimal subtraction scheme. As the study of the diluted Heisenberg mag- 
nets [129] shows, the two-loop approximation is quite precise to describe the 
experimentally observed behaviour of the effective critical exponents. More- 
over, the investigations of effective critical behaviour of the random-site Ising 
model within high-loop order approximations [125] in a semiquantative way 
fully agree with the qualitative picture obtained in the lower approximations 
[117]. Therefore since the effective critical exponents are nonuniversal quanti- 
ties, we consider the two-loop approximation as satisfactory approach for their 
description. 



4-3.1 Isotropic distribution 

The solutions of the system of flow equations (31), with the /3-functions given 
by the resummed expressions (37)- (39) determine the scale-dependent cou- 
plings for the RAM with the random anisotropy axes isotropic distribution. 
The initial values for differential equations usually are chosen in the vicinity 
Gaussian FP describing the background behaviour of system. For the RAM 
these values are restricted by the subspace -u > 0, f > and w < 0. There 
is also an additional condition for the bare couplings: wq/uq = —m. How- 
ever, starting from renormalized coupling values that satisfy this condition, 
the flows fast run into the region with large negative w, where poles appear 
in the integrals for the resummed functions (45). In addition, the existence 
of such a region questions the stability of the Hamiltonian (see Section 3). 
Therefore, we enlarge the region of the initial values of couplings, not neces- 
sarily restricting them to the condition w/u = —m. Since the flow picture for 
the case m = 2 can be analysed without appealing to calculation of princi- 
pal values of integrals in the resummed RG functions and thus is much more 
reliable, we have chosen this case to be presented here. 

In the Fig. 1 we show the flows that start from the initial values with v = 0. 
All flows starting from values of couplings below the separatrix connecting 
FP I and VI run into the region of large negative w, as discussed above (e. g. 
flow 3 in the Fig. 1). In fact they are "runaway" solutions which may indicate 
the possibility of a first order transition. The curves 1 and 2 in the Fig. 1 
are obtained from initial values located above the separatrix connecting FPs I 
and VI. We expect that these flows can not be observed for the RAM because 
they are far away from the condition w/u = —m. The flows starting above 
the separatrix connecting FP I with FP VI are influenced by the unstable 
FP VI and finally attracted by the "polymer" stable FP III. Fig. 2 shows the 
curves of magnetic susceptibility effective critical exponent •jes corresponding 
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Fig. 1. Flows in the space of couplings at w = 0, ?n, = 2. A square denotes the stable 
FP, circles denote unstable FPs and dashed lines show separatrices connecting two 
FPs. Curves 1 and 2 correspond to a flow for the ratio \w/u\ chosen close above 
the separatrix between the FPs I and VI, while the flow 3 starts from a value \w/u\ 
below the separatrix. 



'eff 




Fig. 2. Effective critical exponent 7cflf as a function of the logarithm of the flow pa- 
rameter, m = 2. Curves correspond to flows of Fig. 1 (denoted by the corresponding 
numbers). The dashed lines indicate the values of 7 at the corresponding FP. 

to the flows of the Fig. 1. We do not present results for other effective critical 
exponents, since they demonstrate qualitatively the same behaviour as 7eff. 

Let us now consider the possible scenarios for the effective critical exponents 
if all initial values of couplings are non-zero. Flows starting from physical 
initial values for the RAM are either attracted by the stable FP we rejected 
as unphysical (see the last paragraph of subsection 4.2.1) or have runaway 
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Fig. 3. The effective critical exponent 7eff as a function of the logarithm of the flow 
parameter for isotropic distribution of the local anisotropy axes, m = 2. Dashed 
lines indicate the value of 7 at the corresponding FPs. 

character. However the critical behaviour in the initial phase of the flow can 
be influenced by the physical unstable FPs. At initial values appropriate to 
the RAM a flow can be realized, that is affected by the unstable FP II, which 
describes critical behaviour of pure ?ti- vector model. Thus the effective critical 
exponent corresponding to this flow (Fig. 3) demonstrates regions with the 
critical exponent of the m- vector model (in the given case for vn = 2\ Monte 
Carlo calculations for m = 2 report critical exponents similar to the pure XY 
model [68]. Our theory suggests that the exponents found are effective ones. 



4-3.2 Cubic distribution 

Flows for the cubic distribution of the random anisotropy axes have been 
considered already within the five- loop massive scheme [77], however the be- 
haviour of effective critical exponents has not been studied. They confirm that 
the FPs of the random Ising FP XV and "polymer" FP III models are the 
only stable FPs. Moreover, the flows from initial conditions satisfying u > 0, 
w<0, u + w<0 and v + y > 0, mv + y > according to (21) are never 
attracted by the "polymer" FP III. Therefore in our investigation we restrict 
the initial conditions to this region. Starting with large negative w, the flow 
fast runs into the region, where the resummation method does not work. Thus 
the initial conditions are chosen close to the line w/u = —m. 

As an illustration, let us first consider the trajectories with the initial condi- 
tions M = f = 0. Solutions of the fiow equations in this case are plotted in the 
Fig. 4. The corresponding scenarios of the behaviour of the effective critical 
exponent 7eff are shown in the Fig. 5. The particular features of initial condi- 
tions with M = f = is, that the effective critical exponents are independent 
of the order parameter dimension m since all m-dependent terms in the flow 
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Fig. 4. Flows in space of couplings for u = v = 0. The square denotes the stable 
FP, the circles mean unstable FPs and the dashed lines are separatrices connecting 
two FPs. 



1.35 




-60 -40 

In/ 

Fig. 5. Effective critical exponent j^s ^is a function of the logarithm of the flow 
parameter. Curves correspond to flows of Fig. 4 denoted by corresponding numbers. 
Dashed lines indicate the value of 7 at the corresponding FPs. 

equations drop out. At large values of ratio \w/y\ one can observe flows which 
seem to be not attracted by the stable FP XV at the beginning, but finally 
the flows end in this FP. Therefore effective critical exponents with values 
overshooting the asymptotic ones can be observed (curve 1 in the Fig. 5). If 
flows are strongly affected by the unstable FP V, then before they crossover to 
asymptotic values of the dilute Ising model the critical exponents of the pure 
Ising model (corresponding to FP V) are observed (curve 4 in the Fig. 5). Since 
the effective Hamiltonian (10) with uq = Vq = corresponds to the diluted 
Ising model, in Figs. 4 and 5 the picture obtained for that model [117,125] is 
recovered. 
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Fig. 6. Flows in space of couplings for v = 0. The squares denote stable FPs, the 
circles mean unstable FPs and the dashed lines are separatrices connecting two FPs. 




Fig. 7. Effective critical exponent 7eff as a function of the logarithm of the flow 
parameter. Curves correspond to flows of Fig. 6 denoted by corresponding numbers. 
Dashed lines indicate the value of 7 at the corresponding FPs. 

If we "switch on" the coupling uq (but leave wq = 0) a richer scenario of flows is 
obtained. The typical ones are given in the Fig. 6. Also in this case results are 
independent on the order parameter dimension m, since all m-dependent terms 
drop out of the flow equation. The effective critical exponent 7eff corresponding 
to these flows is shown in the next Fig. 7. There are two FPs V and X that 
correspond to the pure Ising model. So, there exists the possibility that flows 
starting from certain initial values in the space of couplings flrst are affected 
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Fig. 8. Effective critical exponent 7efr as a function of the logarithm of the flow 
parameter for order parameter dimension m = 2. Dashed lines indicate the value of 
7 at the corresponding FPs. 



by the unstable FP V of the pure Ising model, then are attracted by the Ising 
FP X and only then they approach the stable FP XV. It means that within a 
wide range of temperature the critical exponents of pure Ising model may be 
observed before approaching the asymptotic regime (curve 3 in the Fig. 7). A 
characteristic feature of almost all flows is the non-monotonic crossover to the 
asymptotic values. 

The effective critical behaviour, which can be observed if all initial couplings 
are non-zero, depends on m and is described by Figs. 8, 9 and 10 for values 
of spin component m = 2, 3, 4 correspondingly. Apart from the characteristic 
behaviour of effective critical exponent '-/es presented by curves 1-5 of Fig. 5, 
a new scenarios of effective critical behaviour appear. They are presented by 
curves 6 and 7 in the Figs. 8-10, that demonstrate regions with effective critical 
exponent of the pure ?72-vector model (FP II) and the cubic model (FP VIII) 
corresp ondingly. 

In connection with this result, one should mention Monte Carlo simulations of 
the d = 3 RAM with cubic distribution oi the random anisotropy axes [69]. For 
m = 2 and weak anisotropy the results are consistent with a second order phase 
transition with critical exponents, which are close to those of the pure XY 
model. The theory shows that the asymptotic critical behaviour is governed 
by the random Ising critical exponent, FP XV, while critical exponents of pure 
model, FP II, are possible as effective ones in the certain temperature region. 
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Fig. 9. Effective critical exponent 7eff as a function of the logarithm of the flow 
parameter for order parameter dimension m = 3. Dashed lines indicate the value of 
7 at the corresponding FPs. 
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Fig. 10. Effective critical exponent 7eff as a function of the logarithm of the flow 
parameter for order parameter dimension ra = A. Dashed lines indicate the value of 
7 at the corresponding FPs. 

5 Conclusions 



In the present paper we have studied the critical behaviour of three dimen- 
sional random anisotropy magnets. The experimental, numerical and theoret- 
ical investigations performed so far give contradictory results. However, two 
main conclusions can be drawn from the existing data. The random anisotropy 
with isotropic distribution of quenched anisotropy axes destroys the second or- 
der phase transition of the pure system, while for the anisotropic distribution 
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there exists the possibihty of a second order phase transition. 

The standard tool for describing the critical behaviour of different systems, 
the field-theoretical RG, gives results in consistency with these conclusions. In 
particular, the absence of a stable accessible FP for an isotropic distribution 
of the random axes means an absence of a second order phase transition. First 
derived on the basis of a one- loop RG calculation [73], this conclusion found 
its further support in the analysis of the two- [74,75,76] and the five-loop [77] 
massive RG functions refined by resummation. 

The two-loop approximation within the massive RG scheme applied to the 
case, when random axes are distributed along the edges of the tti- dimensional 
hypercube, lead to an answer about the second order phase transition with 
the same set of critical exponents as those of the weakly diluted Ising systems 
[75,76,102]. This was confirmed in a recent five-loop RG calculations [77]. 

Although the majority of the experimental and simulation data show a effec- 
tive critical behaviour, the effective critical behaviour of random anisotropy 
magnets was not studied extensively so far. We used the minimal subtrac- 
tion scheme of the field theoretical RG approach to analyse such a behaviour. 
The calculations were performed within two-loop approximation combined 
with resummation of perturbation theory expansions. The results concerning 
the asymptotic critical behaviour confirm the picture obtained in the massive 
scheme [74,75,76,102,77]. 

Studying nonasymptotic critical behaviour of the d = 3 RAM with isotropic as 
well as with cubic distributions of anisotropy axes we calculated the effective 
critical exponents. The results show that for certain initial conditions the 
effective critical behaviour of the random anisotropy magnets with isotropic 
distribution of random axes can be governed by critical exponent of the m- 
vector model at least for m = 2 in a certain temperature interval. Also we 
investigated effective critical behaviour of systems described by Hamiltonian 
(9) with values of bare couplings differing from those of the RAM. 

In the case of the RAM with cubic distribution of anisotropy axes, various 
scenarios of effective critical behaviour can be observed. The crossover from 
background behaviour described by the Gaussian critical exponents to the 
asymptotic random Ising behaviour demonstrates the existence of temperature 
regions where effective critical exponents coincide with values of the critical 
exponents of pure Ising, cubic and m- vector models. A characteristic feature 
of the majority of the effective critical exponent is the appearance of peaks 
before the crossover to the asymptotic limit. 

In conclusion, we want to attract attention to the importance of determining 
the nonasymptotic critical regime in the investigation of the critical behaviour 
of random anisotropy magnets. As our investigation shows being not close 
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enough to the asymptotic critical regime different values of critical exponents 
can be obtained, they have been shown to be effective ones. 

This work was supported by the Fonds zur Forderung der wissenschaftlichen 
Forschung under Project No. 16574-PHY. Yu. H. and M. D. acknowledge 
French-Ukrainian cooperation Dnipro project and useful comments of Bertrand 
Berche and Taras Yavors'kii. 



A Appendix 



Here, we describe the procedure of the RAM representation in terms of the 
functional integrals. The case of weak anisotropy is of special interest and 
to introduce the controlling parameter D / J an appropriate normalization of 
spins can be used S' = Sv J . Then, Hamiltonian (1) is rewritten in form: 

^ = - E ^^^'r^^'r' - T EI^r^^'r)'- (A.l) 

R,R' '^ "^ R 



It is useful to perform a Fourier transformation for the interaction part of (2) 
and leave the random axis part as it is, after that the Hamiltonian reads: 

H = - E ^{^)S'J'-i. - D'Y.{x^S'^)\ (A.2) 

k R 



where S''k and z/(A;) are Fourier transforms of spin and interaction correspond- 
ingly, D' = D/J. 

For the fixed configuration of the local random anisotropy axes {x-r} the 
configuration-dependent partition function is written 

Z({£r}) = Sp e-'^^, (A.3) 



where the Hamiltonian Ti is given by (A.2), j3 is the inverse temperature and 
the trace means integration over the surface of TTi-dimensional hypersphere 
with unit radius: Sp(. . .) = J dS'ni- ■ •)^(I'S''r| — !)• 

To avoid taking trace of products of spins one can apply the Stratonovich- 
Hubbard transformation to expression (A.3) and arrive at the following form 
of the partition function: 



Z{{xr}) ~ Sp y {d(f) exp { - - ^ ^^Jm ''^'''^ ^ 
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E 'S'Wr + PD'Y^ixnS'uf}, (A.4) 

R R 

where 0k is the Fourier transform of the m-component field variable 0r,, in- 
troduced by the Stratonovich-Hubbard transformation. The integral in (A.4) 
means a functional integration in the space of the field variables 0. Here, we 
omit the normalization coefficient before the integral since it does not alter 
the behaviour of system near critical point. 

Now the trace operation for the spins concerns only the last two terms of 
the exponent in (A.4). Expanding the exponential function in a Taylor series, 
performing the integration with use of generalized m-dimensional spherical 
coordinates and finally reexponentiating the obtained series, one gets: 



Zi{xn})-J{d$)exp{-^Y. 



kP Er |0Rp I -, 



, f3u{k) 2m 

ER|0R|'(^ + 4 + /3D'(m + 6)) , 



4m2(m + 2)(m + 4) 



2 




m? (m + 2) (m + A) J' 

Here, we omit all terms independent on since they give only a shift of the 
free energy. Expression (A. 5) has been obtained under assumption of a weak 
anisotropy D' << 1, thus the only terms in first order in D' were taken into 
account. Furthermore, there is no interest in terms of order higher then fourth 
order in 0, since they are irrelevant in the RG calculations of the critical 
properties of the model. They may be relevant when multicritical behaviour 
is studied. 

Since the small wave lengths are relevant for the critical behaviour, it is 
possible to use the expansion of the interaction Fourier image at small k: 
{j3v{k))^ = Co + &ok^ + . . .. Then transforming back to real space and pass- 
ing to the continuous limit by replacing summation over R by the familiar 
(i-dimensional integral one obtains: 



Z({£r}~ /(rf0)e^(^^'<^^), (A.6) 

where 



n{xj,Jj,) = -Jd''R{^ 



roi0Rr+ 



V0r' 



-Di (0r£r +t;i|0R| + 



;^|0Rr (^r^^r) +•••}• (A.7) 
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Here, the field was normalized to get coefficient for gradient term equal to 

unity ana ro - ir Cq - -1,1- „r„ , o^ ^ ' ^l " hn^m2fm+2Vm+4^ ' ^l " ft„^mfm- 



bo V '^ m^ m(m+2)^/' -^ bo m2(m+2)(m+4) ' ^ 6o m(m+2) ' 

6o^m2(m+2)(m+4) " 
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